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Abstract. Algebraic space-time coding allows for reliable data exchange 
across fading multiple-input multiple-output channels. A powerful tech¬ 
nique for decoding space-time codes is Maximum-Likelihood (ML) de¬ 
coding, but well-performing and widely-used codes such as the Golden 
code often suffer from high ML-decoding complexity. In this article, a 
recursive algorithm for decoding general algebraic space-time codes of 
arbitrary dimension is proposed, which reduces the worst-case decoding 
complexity from 0(|5|” ) to 0(|5|"'). 


1 Introduction 

Algebraic Space-Time (ST) coding is a powerful technique for reliable data ex¬ 
change across fading Multiple-Input Multiple-Output (MIMO) channels, which 
makes use of available multiple antennas at the transmitter and receiver, as well 
as multiple channel uses to introduce redundancy in the transmitted information. 
Data transmission across a MIMO channel can be modeled as 

knrXT — XntXT + Nn^xT, (1) 

where the subscripts rit, Ur and T denote the number of antennas at the trans¬ 
mitter, receiver, and the coding delay, respectively, Y and X are the received 
and sent matrices, iL is a random complex matrix modeling Rayleigh fading, and 
TV is a noise matrix whose entries are complex Gaussian with zero mean. In this 
paper we restrict ourselves to the fully symmetric case where rit = Ur = T = n 
for some n, and thus omit the above subscripts. 

First introduced two decades ago, ST codes have been the subject of exten¬ 
sive study, and various algebraic criteria have been derived for ensuring reliable 
performance. The complexity of the decoding algorithm is a major issue in prac¬ 
tical implementation. In digital video broadcasting [I] , high decoding complexity 
prevents the use of high-rate space-time codes, since the complexity of linear de¬ 
coders such as the sphere decoder grows exponentially in rank. 
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For the remainder of this article, we let vec : Mat(TO x n, C) —>■ C"^" be the 
vectorization transformation and we denote by S the signaling alphabet {e.g. 4- 
QAM) used. Write vec{HX) = Gs, where s G 5^. Decoding amounts to solving 

argmin ||F — HX\\^p argmin ||vec(F) — GsH^, (2) 

XeX 865 '= 

where || • ||f and || • || denote the Frobenius and Euclidean norm, respectively. 
The latter search can be performed with the help of a sphere-decoder [5] , whose 
decoding complexity is upper-bounded by 

Cleverly constructed algebraic ST codes known as fast-decodable codes reduce 
the complexity of ML-decoding by parallelizing the ML search (see e.g. E])- 
However, ST codes exist that, although they do not allow for a reduced ML- 
decoding complexity, exhibit other desirable properties. An example of such a 
code is the Golden code [1], incorporated in the IEEE 802.16 (WiMAX) standard. 
The main algorithm of [5] reduces the decoding complexity of the Golden code 
from 0(|5|^), corresponding to the complexity of ML-decoding, to 0(|5p), while 
maintaining nearly-ML performance. The algorithm is specific to the Golden 
code, but has been generalized to the 3x3 and 4x4 perfect codes in, respectively, 
[5] and [7]. The main contribution of this article is a recursive algorithm for 
decoding arbitrary n x n algebraic ST codes generalizing those of mm, which 
reduces the worst-case decoding complexity of 0(151" ) down to 0(|5|"). 

We review how algebraic ST codes are constructed from central simple al¬ 
gebras in Section [5] Our algorithm for decoding general n x n algebraic ST 
codes is illustrated in Section [3l wherein we also discuss its complexity. We then 
briefly consider the case n = 2 and the Golden code as a special example in Sec¬ 
tion H) In Section [5l we study the 4x4 perfect code in more detail and provide 
an elementary statistical analysis of the involved quantities. While preliminary 
simulations seem to indicate that our algorithm has performance which is close 
to ML-decoding, time and space constraints prevent us from providing proper 
error rate simulations in this current draft. Detailed error rate simulation results 
will appear in a subsequent version of the paper. 


2 Space—Time Codes from Central Simple Algebras 

Let LjK be a Galois number field extension of degree n, and let A be a K- 
central simple L-algebra of rank n over L, which we treat as a right L-vector 
space. For the rest of the article, we fix compatible embeddings of K and L into 
C, and identify K and L with their images under these embeddings. To construct 
matrices to serve as codewords, we consider an injective algebra homomorphism 

p : A —?► Mat(n, C). (3) 


Definition 1. An algebraic space-time code is a finite subset X of p{A) or its 
transpose p{AY. 


As a general reference for applications of central simple algebras to space-time 
coding, we recommend [S], and as a catch-all mathematical reference for central 
simple algebras, we will use [9l Chapter IV]. There are many properties that 
can be desirable for an algebraic ST code, such as full-rank codewords, non¬ 
vanishing determinants, fast-decodability, balanced energy across transmitters, 
and optimal diversity-multiplexing tradeoff, which can be achieved by demanding 
some additional algebraic requirements of the involved structures. 

Let us fix a basis {ei,..., e„} of .4 as a right L-vector space. As an example 
of the above, consider the left regular representation (LRR), defined as follows. 
For fixed x £ A, the map given by left multiplication y xy for aX\ y £ A 
is right L-linear and compatible with algebra multiplication, and thus allows 
us to represent x as a complex matrix p{x) in the above basis. For notational 
convenience, we will below mostly consider the transpose of this representation. 

To give a concrete example, consider the case where L/K is a cyclic extension, 
whose Galois group Gal(L/Ar) = (ct) is generated by a single element. Fix some 
7 S AT^ and consider the central simple algebra 


n—1 




( 4 ) 


where the multiplication is determined by 


Xcj = ej(j^{\) for all A G L, and 


( 5 ) 

( 6 ) 



If p denotes the LRR, the explicit structure of the codewords of X C p{Ay 
can be found in, for example, [8]. More generally one can consider the left 
regular representation of the crossed-product algebra A((p) for any 2-cocycle 
if : Gal(L/iF)2 ^ 0 Ghapter IV]. 

Let us return to the general case oiL/K a. Galois extension of degree n, and A 
any AT-central simple L-algebra of rank n over L. Gonsider any injective algebra 
homomorphism p : A ^ Mat(n,C). The following decomposition of the matrix 
p{xy for X £ A and its vectorization are central to our decoding algorithm. 

Given x G M, we can represent it as x = ^ 

{uji,... lOJn} to be a basis for L as a AT-vector space, then we can write x as 



where Xji £ K (7) 


Passing to the transpose of the representation p gives 



( 8 ) 


where we have used the fact that pixji) is a scalar matrix and thus commutes 
with everything. If we define Xi = Pi^jY then vectorizing p{xY 

gives, by basic properties of the vectorization operator, 


n n 

vecipixY) = ^{In ® vec(Xj) = ^(/„ 0 p{uJiY)r,Si 
2=1 2=1 


(9) 


where Si = [xu ■ ■ ■ Xnif and Fi G Mat(n^ x n, C) is the unique matrix satisfying 
Fi = [r'fj • • • Ff^^Y where FjiSi is the column of Xi. Note that the above also 
provides us with a decomposition of the generator matrix G of the code, as the 
block matrix 

G = [{In 0 p{uJlY)Fi • • • (/n 0 p{uJnY)rn] (10) 

It is worth pointing out that in the case of the cyclic algebra constructed in the 
above example, the matrices Fi are the same for all i. 

The elements Xji G K are the information symbols we wish to decode, and 
one should think of them as belonging to our signaling alphabet S. To introduce 
some discrete structure, one usually assumes that xji G Ok, where Ok denotes 
the ring of integers of K. The basis wi,..., a;„ is often then chosen to be a basis 
for Ol as an Oi^-module (assuming this is possible), or more generally, a basis 
of an ideal I C as an O^-module. In the case where I = (a) C Ol is a 
principal ideal, the effect of considering this ideal amounts to replacing uJi with 
ujia in the above expressions. Such an a is often called a shaping element. 


3 The General n x n Case 

In this section we present the main result of the paper, which is a recursive 
decoding algorithm for n x n algebraic ST codes with complexity 0(|5|"). Infor¬ 
mally, the algorithm uses the decomposition of vec(X) of Section [5] to split up 
the information symbols contained in an n x n codeword into n vectors, each 
containing n symbols. Based on how well-behaved certain channel sub-matrices 
are, we then estimate one of these vectors, requiring a search over the space 
|iS|". This estimate is subtracted from the received signal and we then proceed 
by recursion. 

Our algorithm performs n searches over the space |5|", rather than the single 
search over the space |iS|” required by the naive brute-force algorithm. Hence 
it is clear that the complexity is OdiSI”), when measured in terms of the size 
of the space we must search over when solve the below “arg min”-type problem. 
However, one pays for the reduction in complexity in this operation by com¬ 
puting several covariance matrices, evaluating their determinants and condition 
numbers, and applying a lattice quantizer. Hence, as with the main algorithm 
of [3 which we are generalizing, there is a non-trivial amount of preprocessing 
to be done. This algorithm also generalizes the main algorithm of [7], wherein 
the authors break up the 4x4 perfect code into two groups of eight symbols, as 
well as the approach in which studies the 3x3 perfect code. 



Let LfK, A, and p be now as in Section O and consider an algebraic ST 
code of the form X C p{Ay . Given the channel equation ([T|) for a channel 

finite 

matrix H = {hki)i<k,i<n, we first vectorize the expression to get vec(y) = 
vec{HX) + vec{N). By similar arguments as in Section[2l we have 

n n 

yec{HX) = ^ O {Hp{uJiy))riSi = ^ HiSi, (11) 

i=l i=l 


where, using the same notation as in Section [5J Si = [xu ■ ■ ■ is the vector 

containing n of the symbols we wish to decode, and Hi = {In®{H p{u)iy))ri S 

n 

Mat(n^xn, C). Setting y := vec(V), we can rewrite equation ([T]) as y = ^ HiSi+ 


vec(X). Decoding the vector y now amounts to solving 


arg min 
(si,...,s„)e(S") 


-ZH,. 


( 12 ) 


The idea of the decoding algorithm is to first decode some , use this infor¬ 
mation to recursively decode some Si^, and so on until ultimately decoding Si^ 


in the step. Define Gi = 


. Let k(-) denote the condition 


Hi-- - H,-- - Hr, 

number of any matrix. The algorithm chooses which vector of symbols to decode 
first based on the following decision. 


Decode Si^ first, where ii = arg max 


det(GlG,) ] 

<GlGy J 


(13) 


As is noted in [^, the determinant of the above covariance matrix measures 
the instantaneous SNR of the corresponding linear system and thus should be 
large, and the condition number measures the accuracy zero-forcing approxima¬ 
tions and thus should be small. To decode Si^ we will express the other symbols as 
a function of Sir y for this function to be as tolerant to interference and noise 
as possible, one should maximize the above determinant and minimize the above 
condition number. Thus we suggest maximizing their ratio as a compromise. 

Reindex the terms HiSi such that = 1 and define ti = [s^ • • • Then, 
the expression involved in equation (USD can be rewritten as 


n 

y-'^ HiSi 

2 

2 = 1 



y- HiSi - y^ HjSj 


2=2 


= {y- iLiSi )^(/„2 - Gi{G\Gi)-^G\){y - Hisi) 

+ (D('Si) — ti)^G|Gi(fi(si) — ti), 


(14) 

(15) 

(16) 


















(17) 


where Gi = [iJ 2 • ■ • -ffn], and 

fi(si) = {G\Gi)-^Gi{y - Hisi). 

Here we are defining ti(si) to be an approximation of ti given si. The ap¬ 
proximation is essentially chosen to minimize the value of the quadratic form 
u^G|Giu, which in turn should minimize the expression in (na. The success of 
this process will depend on the determinant and condition number of G|Gi in 
the manner stated above. 

Let 13 be a quantizer for the signaling alphabet S. Since ti(si) is not neces¬ 
sarily a constellation point, we now decode an estimate 

fi(si) = i3(ti(si)), (18) 

and substitute this function into (na, which now only depends on si and reads 

argmin y — TLiSi — Giti(si) . (19) 

Si 

We can thus compute an estimate si of si. Having estimated 5 i, we can go back 
to (na and replace y with y — HiSi. Using the same procedure and reindexing, 
we can solve for S 2 - The algorithm is done after n steps. We present the algorithm 
in pseudocode below. 


4 The 2x2 Case 


As a special case of the above scenario, set n = 2 and let L/K be a degree 2 
real extension of number fields, that is L = K{'/d) for a squarefree d € Z>i, 
and G\k = 1. Let Ok and Ol = Gic[w] be the ring of integers of K and L, 
respectively. Necessarily, the Galois group Gd\{L/K) = (a) oi L/K is cyclic of 
order 2 , where a is dehned by a : Vd i—> —Vd. 

Choose 7 € such that 7 ^ Nmj^/^(L^) and I 7 P = 1, and define the 
quaternion algebra 

Q = (d, 7 )/f = L © eL, ( 22 ) 

where = 7 . For a shaping element a G , the ST code constructed from the 
transpose of the left regular representation of Q is a subset 


A C 

finite 


r(a) 


Xi + X2UJ X3 + X4IJJ 

7(3:3 + a:4CT(a;)) xi + X 2 (j{uj) 


Xi G Ok, 1 < * < 4 


(23) 


Given the transmission model in o for a channel matrix H = [ ], 

vectorizing the received matrix results in the equivalent equation 


vec(F) = vec{HX) + vec{N) = Hisi + H 2 S 2 + vec{N), 


(24) 


where si = [xi X2Y', S2 = 3:4]*, and 







Algorithm 1 Recursive Decoding of ST Codes 


1: Input: Y, H 
2: Output: si,..., Sn 

3: Step 1: 

4: for i = 1,..., n do 

5: Compute Ri = (/„ (g> (Rp(cJi)‘))ri 

6: Define T-L = {Hi,..., H„}, T = {I,... ,n\ 
7: Compute y = vec{Y) 


8 : 

9: 

10 : 

11 : 

12 : 


Step 2: 

while H 7 ^ 0 do 
for i do 

Define Gi = [Hi ■ ■ ■ Hi ■ ■ ■ R„] where Hj € H, j ^ i. 


Compute 


ii = arg max 
i€I 


det(GlG0 ] 

^{GlGi) J 


13: Ypdate H = H\ {Hi,}, I = T\{ii} 


14: Step 3: 

15: Compute 

Sii = arg min \y - Hi,Si, - Gi,Q{{Gl Gi,)~^Gi,{y - 

Si,es^ 11 

16: Update j/= 2 /— RiiSii 

17: goto Step 2 


Hi, Si,)) 


2 


( 20 ) 


( 21 ) 


cxhii 

a{a)'yhi2 


arhii 

a{aT)'-fhi2 

a{a)hi2 

cxhii 

, H2 = 

a{aT)hi2 

arhii 

ah2i 

a{a)jh22 

aTh2i 

(t(q!t)7/122 

a{a)h22 

ah2i 


(j{aT)h22 

aT/i2i 


The relevant covariance matrices are 


ijjiJi = 


and h\H 2 = 


a^j NmL/K{uj)b 

(NiriL/Kii^)b)* ttui 


where the superscript * denotes complex conjugation, and 

a= ||aij(i)||2 + ||a(a)i7(2)||2^ 

+ \\a{auj)H^'^''\\^, 

b = (aH‘'^\aia)HG)) + (a(a)-iH^‘^\aHG)). 


(25) 


(26) 


(27) 

(28) 
(29) 


where is the column of H. These covariance matrices exhibit a very 
simple form due to the assumptions that L is a real extension of K and jyp = 1. 
The determinants and condition numbers can easily be computed to be 








Hi Hi 

hIH2 

det(-) 

a 2 -| 6|2 

al - 1 Nm4,/;f(a;)6p 

k(-) 

a-|-|b| 

a^ + \ Nmj:,/K(w)h| 

“-FI 



Thus, for example, if 

^ det(g|gi) 

k{hIH 2 ) ~ k{hIHi) 

then we define the following function of s £ 5^ 

his) = iHlH2)-^H2iy - His) (31) 

and subsequently estimate Si and S 2 by computing 


Si = argmin||y - 77is - i 72 S 2 (s)|p, S 2 = argmin ||?/- i/isi - 772 s|| (32) 

sGS2 sGS2 

For the Golden code [3], which is the code of interest in [S], the quaternion 
algebra used for code construction is Qg = (5,i)Q(i) with uj = ■ In this 

special case it is easy to see that the covariance matrix with the larger determi¬ 
nant also has the smaller condition number. Thus one can reduce to the decision 
metric to simply computing argmax^^j^ 2 {det(i?/i?i)}. 


5 The 4x4 Perfect Code 

In this section we illustrate the decomposition of Section [5] for the 4x4 perfect 
code. Let L/K = Qii,0)/Qii), where 0 = Ci 5 + Ci 5 ^ = 2 cos (27r/15). L/K is 
a cyclic extension of degree 4, and its Galois group is Gal(L/iF) = (tr), where 

3 

O' ■ Ci5 + Ci 5 + Gonsider the central simple algebra .4=0 ejL as 

3=0 

defined in Section[5] If cc £ .4 and p is the left regular representation, a codeword 
is of the form X = p(a:)* and the corresponding decomposition of vec{HX) is 

4 

vec{HX) = ^(/4 ® iHpiu3i)))rs, (33) 

where the oji are given explicitly by wi = (1 — 2>i) -\- a ;2 = (1 — 3z)d -I- 

ujz = —i-b (—3-|-4i)0-|- (1 — i)9^, and a ;4 = (-1-1-*) — W-\-+ 9^ and the matrix 
T = [r 4 •.. rlf is given by 


4 0 0 O' 


'0 10 0' 


'0 0 10' 


'0 0 0 1' 

C 

O 

o 

o 

, r 2 = 

10 0 0 

, = 

0100 

, A = 

0 0 10 

0 0 7 0 

0 0 0 7 

1000 

0 10 0 

0 7 0 0 


0070 


0 0 0 7 


10 0 0 


( 34 ) 
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Fig. 1. The empirical distribution of mini{Ai(G'tG'i)} during the first round of recursive 
decoding. 

In Fig. 1 we plot the empirical distribution of mini{K(G|G'i)} averaged over 
10® channel matrices. As in [2, we notice that the minimal condition number is 
quite well-behaved, hence it is highly likely that one of the four linear systems in 
question is well-conditioned and thus we can expect decoding to be successful for 
a high percentage of channels. In contrast with [3, it is no longer true that the 
index minimizing the condition number is the one maximizing the determinant. 
While we are not deciding based solely on the condition number, the index 
minimizing the condition number is highly correlated with the one maximizing 
the ratio, having correlation coefficient ~ 0.6470. 



Fig. 2. The empirical distributions of maxi{det(GlGi)} and 
maxi{det(G|Gi)/«:(G|Gi)}, respectively, during the first round of recursive de¬ 
coding. 


In Fig. 2 we compare the empirical distributions of maxi{det(G|Gi)} and 
maxi{det(G|Gi)/K(G|Gi)}. The similarity of the histograms is a symptom of 
the high correlation coefficient of ~ 0.9095 between the i maximizing the de- 


















terminant, and that which maximizes the ratio. Thus the decision metric in the 
algorithm essentially amounts to selecting the index maximizing the determi¬ 
nant, with the condition number selecting the better-conditioned linear system 
among several with similarly-sized determinants. 

6 Conclusions and Future Work 

In this article, we have proposed a recursive decoding method for algebraic space- 
time codes of arbitrary size, which has complexity 0(|iS|”) compared to the 
worst-case complexity of 0(|5|" ) of ML-decoding. The algorithm is based on a 
decomposition of the representation of the corresponding central simple algebra 
into a sum of n matrices, each of which encodes n information symbols. These 
groups of symbols are then decoded recursively, based on certain properties of 
equivalent channel sub-matrices. As basic examples, we have studied the relevant 
algebraic decompositions for 2 x 2 codes and the 4x4 perfect code, as well as 
provided a first statistical understanding of the algorithm for the latter. 

The most immediate task at hand is that of experimental error rates for 
various algebraic space-time codes, comparing our algorithm with ML-decoding 
and with the algorithms proposed in m- This simulation results are forth¬ 
coming and will appear in a later version of this paper. One promising avenue 
of future work is comparing and possibly combining this approach with that of 
fast-decodability. The two approaches appear mutually exclusive, and thus could 
be combined to further reduce decoding complexity. 
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